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Abstract 

This paper presents the simulation of the GLAST high energy gamma-ray 
telescope. The simulation package, written in C-I--I-, is based on the Geant4 
toolkit, and it is integrated into a general framework used to process events. 
A detailed simulation of the electronic signals inside Silicon detectors has 
been provided and it is used for the particle tracking, which is handled by 
a dedicated software. A unique repository for the geometrical description of 
the detector has been realized using the XML language and a G-I-+ library 
to access this information has been designed and implemented. A new event 
display based on the HepRep protocol was implemented. The full simulation 
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was used to simulate a full week of GLAST high energy gamma-ray observa¬ 
tions. This paper outlines the contribution developed by the Italian GLAST 
software group. 


1 Introduction 

The Gamma-ray Large Area Space Telescope (GLAST) is an international mission 
that will study the high-energy phenomena in gamma-rays universe GLAST 
is scheduled for launch in mid 2007. 

GLAST is instrumented with a hodoscope of Silicon planes with slabs of con¬ 
verter, followed by a calorimeter; the hodoscope is surrounded by an anticoinci¬ 
dence (ACD). This instrument, called the Large Area Telescope LAT, is sensitive 
to gamma rays in the energy range between 30 MeV and 300 GeV. The energy 
range, the field of view and the angular resolution of the GLAST LAT are vastly 
improved in comparison with those of its predecessor EGRET (operating in 1991- 
2000), so that the LAT will provide a factor of 30 or more advance in sensitivity. 
This improvement should enable the detection of several thousands of new high- 
energy sources and allow the study of gamma-ray bursts and other transients, the 
resolution of the gamma-ray sky and diffuse emission, the search for evidence of 
dark matter and the detection of AGNs, pulsars and SNRs. A detailed description 
of the scientific goals of GLAST mission and an introduction to the experiment 
can be found in |5]. 

GLAST is a complex system, and detailed computer simulations are required 
to design the instrument, to construct the response function and to predict the 
background in the orbit. To accomplish these tasks an object-oriented C-I--I- ap¬ 
plication called Gleam (GLAST LAT Event Analysis Machine) was adopted and 
implemented by the GLAST LAT collaboration. A brief description of Gleam 
could be found in 0. Its structure is described in figure H An entire week of the 
gamma-ray observations by GLAST-LAT was simulated using Gleam in order to 
develop and test GLAST LAT scientific analysis softaware. 


2 Simulation and Reconstruction Software 

The GLAST off-line software is based manly on Gaudi, a C-I--I- framework, orig¬ 
inally developed at CERN Gaudi manages the loop of particles to be sim¬ 
ulated, then a series of algorithms are applied to each of them to get the result 
of the complete simulation and reconstruction chain. The Source Generation is 
the first algorithm called within the particle loop. Its task is to generate particles 
according to certain characteristcs. This algorithm must store the information on 
the temporal and spectral behaviour of the source, as well as on the orbital char¬ 
acteristics of GLAST. It provides a user interface to produce additional incoming 
particles and is responsible for setting the current time, the particle energy, direc¬ 
tion, and type. Within this package a series of default sources are implemented. 
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GLEAM Application 
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Figure 1. General scheme for simulation and reconstruction within the GLAST off-line software 
framework 


They include source for testing purposes as well as the description of astrophysical 
spectra and the expected particle and albedo gamma backgrounds. An extension 
of this framework has been implemented for simulating transient sources such as 
Gamma-Ray Bursts (GRB). It can be used for studying the capability of GLAST 
for the observation of rapid transient fluxes in generally. 

The algorithm which is responsible for generating the interactions of particles 
with the detector is based on the Geant4 MonteCarlo toolkit which is an Object 
Oriented (00) simulator of the passage of particles through matter. Its applica¬ 
tion areas include high energy physics and nuclear experiments, medical science, 
accelerator and space physics. Within the Gleam application the simulation is 
managed by the Gaudi algorithm G4Generator. 

Figure [21 shows an event generated using Geant4 within the GLAST LAT ex¬ 
periment. 

The next algorithm to be applied is the digitization which transforms the hits 
generated by the Monte Carlo simulations into the signal as read by the electronics. 
To implement a detailed digitization of the Tracker system a full simulation code 
has been developedjj). Starting from Monte Carlo hits in the Silicon detectors, the 
current signals induced on each strip are evaluated and are converted into voltage 
signals using the transfer function associated to the detector electronics, taking 
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Figure 2. High energy gamma-ray interacting with the GLAST LAT detector 


into account the detector noise as well as the noise associated to the electronics. 
The fired strips and the time over threshold are then determined. 

The signals in the Tracker are then analysed by the reconstruction package. It 
generates a series of clusters, that are used to find and fit the best track candidates. 
This last procedure is done using alternative pattern recognition algorithms and 
a Kalman Filter based algorithm. Finally, using the best track found, another 
algorithm finds the best vertex candidate for gamma events. 

Although it is not part of the simulation, the visualization package is essential 
for the use of the simulation itself. A new version iHiini of the event display based 
on the HepRep protocol was developed and integrated in the offline software. The 
figure 121 shows a recent FRED-based event display of GLAST. 


3 Simulating the GLAST LAT High Energy Sky 

Using Gleam it’s possible to compute and store time, direction and energy of 
each incoming simulated and reconstructed gamma ray. The Data Challenge One 
(DCl), organized by the LAT Collaboration from December 2003 to February 
2004, represented the first opportunity to test the complete simulation chain, and 
the first attempt to perform scientific analyses on simulated data. For DCl, only 
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Figure 3. GLAST LAT event display based on FRED 


the gamma-ray sky was simulated, while the cosmic-ray flux (about 10^ times 
greater) was modeled separately for development of the background rejection al¬ 
gorithms. These algorithms were then applied to the simulated gamma-ray data. 
The gamma-ray sky simulated for DCl was quite rich: a variety of sources was 
included. The simulation software takes into account the relative fluxes and com¬ 
putes from which source the next photon arrives. A map of the simulated sky is 
reported in fig. 01 



Figure 4- Full simulated gamma-ray sky seen by GLAST LAT 


Based on the simulated data, several groups tried to develop the best algorithms 
to detect unknown gamma-ray sources. GLAST observations contain photons 
from astrophysical sources, convoluted with the spatial and spectral instrument 
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response. Moreover in most astronomical gamma-ray images a large fraction of 
sources is near the detection limit. Thus a careful statistical treatment is needed to 
determine their existence and properties. Many tools (parametric methods) need 
a priori model to fit the data and estimate their parameters. No model or hypote- 
sis on the data are requested by the wavelet method Through an iterative 
procedure the wavelet method allows a blind search of gamma-ray point sources. 

Another goal for the DCl was to test the possibility to trigger Gamma-Ray 
Bursts using only LAT information. Several groups within the LAT collaboration 
are prototyping trigger and alert algorithms for detecting transient signals in DCl 
data, and different algorithms were studied^J. Different algorithms were success¬ 
fully applied for searching for transient signals in DCl data. Bright bursts (with 
fluence greater than 10“® erg/cm^ between 20 keV and 1 MeV) can be detected 
with simple algorithms. Further studies will include the particle background, and 
the possibility to implement an on-board LAT alert algorithm. All of these items 
will be addressed for the next Data Challenge (DC2), in which one month of 
simulated data will be produced. 


4 Conclusions 

The Gleam simulation program has been developed in the last few years and now 
it’s ready for simulating the full GLAST satellite and is being used for deriving 
the final instrumental parameters and for generating a full set of events for the 
developing of scientific analysis software. 
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